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ABSTRACT 



Context. The discovery of new massive star clusters and massive stellar populations in previously known clusters in our Galaxy by 
means of infrared studies has changed our view of the Milky Way from an inactive to an active star-forming machine. Within this 
scenario, we present a near-infrared spectrophotometric study of the stellar content of the compact H II region Sh2-152. 
Aims. We aim to determine the distance, extinction, age, and mass of Sh2-152, using for the first time near-infrared stellar classification 
for several sources in the region. 

Methods. Using our near-infrared (J, H, and K s ) photometry and the colour-magnitude diagram for the cluster field, we selected 13 
bright stars, candidate members of the reddened cluster's main sequence, for H- and K- spectroscopy and spectral classification. This 
near-infrared information was complemented with an optical spectrum of the ionizing central star to confirm its spectral nature. 
Results. From the 13 spectroscopically observed stars, 5 were classified as B-type, 3 as G-type, 2 were young stellar objects (YSOs), 
and 3 remained unclassified (because of the poor data quality). The cluster's extinction varies from A Ks = 0.5 to 2.6 magnitudes 
(Ay = 4.5 to 24 magnitudes) and the distance is estimated to be 3.21 + 0.21 kpc. The age of the cluster is younger than 9.4 Myr and 
the lower limit to the total mass of the cluster is (2.45 ± 0.79) • 10 3 M G . We compare the number of ionizing photons emitted from the 
OB-type stars with the Lyman continuum photons derived from the radio observations and conclude that both quantities are consistent 
for the central region of Sh2-152. In contrast, the main ionizing source of the lower region remains unidentified. 



Key words. Infrared: stars 
spectroscopic. 



1. Introduction 



Stars: early-types, Hertzsprung-Russell and C-M diagrams, mass function - Techniques: photometric, 



Thank s to infrared s tudies such as D ENIS ( lEpchtein et al.L 



1 19971). 2MASS dSkrutskie et all I2006I). GLIMPSE 



dBeniamin etail l2003h . and UKIDSS dLawrence et all l2007l) . 
we have been able to observe regions of high visual extinction, 
previously impossible to achieve with optical instruments, and 
discover tha t our Galaxy is an active star-formation machine 
dKgeEHool), fed by massive clusters and their massive stars. 

Some of these massive obj e cts, fo r example the Arches- 
Quintuplet cluster dFiger et all I1999I). RSGC1 dFiger et all 
l2006llDavies et alll2008l). RSGC2 dPavies et alll2fj07l) . RSGC3 
dClark et all l2009l I Alexander et all 



2009I) and Alicante 8 (a.k.a. 



RSGC4, iNegueruela et al.ll2010l) are ne wly discovered mas sive 
clusters; but others such as Cy gOB2 dKnodlseden 1 l2000h or 
Westerlundl dClark et all 120051) . are previously known clus- 
ters with a newly detected massive stellar population. This 
increase in massive cluster discoveries has been comple- 
mented wi t h new cat alogues of cluster c a ndidates, for exampl e 
iBica et all d2003allbl) iDutra et al.l d2003l) . iMercer et al.l d2005l) . 
and Froebrich et al r d2007l) . If we also consider that hundreds of 
Galac tic massive clusters remain unknown dHanson & PopescuL 
120081) . it is easy to grasp the need to generate projects dedicated 
to finding and characterizing these new massive objects. 



This re-evaluation of our ideas of the Milky Way motivated 
us to develop the MASGOMAS (MAssive Stars in Galactic 
Obscured MAssive clusters) project. Within this project, we 
have completed a photometric near-infrared broad band (J, H 
and Ks ) catalogu e of 44 cluster candidates extracted fro m pre- 
vious catalogues dDutra & BicaL 1200 ll IBica et a n i20033bh and 
we have continued with a medium resolution multi-object spec- 
troscopic follow-up for 9 of the most promising candidates. A 
detailed descr iption of MASGOMAS is given in the first paper 
of the project dMarfn-Franch et alll2009l) . 

One of our candidates is the star-forming compact HII re- 
gion Sh2-152 (a 2 ooo = 22 h 58 m 45 s and 5 2 ooo = +58°46'50"). 
This region is located in the Perseus arm, close to the Galactic 
plane (/ = 108!76,fo = -0°95) and stands out from the back- 
ground in the Spitzer 3.6 fim image shown in Fig. Q] Sh2- 
152 contains two IRAS sources: IRA S 22566+5830 and IRAS 
22566+5828 dKleinmann et allll986l) . The first is located in the 
central zone of the clus ter and woul d be excited by an Q 9 V star 
dCrampton et all 1 19781) . labelled bv lRusseil et alj d2007l) as star 
number 4, and the second IRAS region, located ~ 2 arcmin SE 
from the central section, appears more extended and less popu- 
lated than the central IRAS 22566+5830. 

In the Sh2-152 field of view we find five masers, all of 
them presented in the false-colour LIRIS image (see Fig. 0: 
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one 6.7 GHz methanol ( Szvmczak et all l2000t) . three water 
(iPalad et all I1993L lHariu et all Il998h .~and one hydroxyl (OH) 
dWouterloot & Habingl 1 19850 . These masers are indicative of 
star formation a nd in the case of the hydro xyl one, of high-mass 
star formation dZinnecker & Yorkd [2007b ." Its presence makes 
this cluster an interesting candidate as a sti ll young and star- 
forming active region. Following this idea, IChen et aTl (120091) 
used K', and H2 photometry to detect H2 emission in Sh2-152, 
implying that the cluster hosts ongoing star formation. 

However, even if this region has been the object of near- 
infrared observations its stellar content has not yet been 
characterized spectroscopically. Only the ionizing central star 
has been optically observed and spectro scopically classified 
dCrampton et allll978llRusseil et alll2007l) . 

Anoth er interesting po i nt abo ut Sh2-152 is its distance and 
extension. ICrampton et all d 19781) . reported 3.6 kpc as the first 
estimate of the cluster's distance using the spectroscopic paral- 
lax to the ionizing star. Later, in a series of studies observing 
H2O masers, NH3 emission lines, and both l2 CO(J = 2 — > 1) 
and n CO(J = 3 -» 2) transitions, IWouterloot & Walmslevl 
d!986l) estimated the distance to Sh2-152. Assuming that Sh2- 
152 is located in the Perseus arm, between 2 and 6 kpc or an 
average distance of 3500 pc, they confirm this conjecture first 
with the H2O maser observa tion, and then, with the ammonia 
core dWouterloot et all 119881) . In any case, in the third work of 
this series, IWouterloot et al.l (fl99l . determined a (kinematic) 
distance of 5.3 kpc through radial velocity measurements of 
12 CO(7 = 2 -> 1) and 12 CO( 7 = 3^2) lines for IRAS 
22566+5830. In the same year, lHariu et al.l d 1993b presented 
consistent results, based on observations of ammonia clumps, 
with the initial assumption of a distance of 3.5 kpc for Sh2-152. 

On the basis of lon g-slit optical s pectro scopy for the cen- 
tral region of Sh2-152. iRusseil et al.l d2007l) then calculated a 
distance of 2.39 kpc. This distance was derived from the spec- 
tral classification of the ionizing central star (the aforementioned 
star number 4), although we note that the position of the star is 
inaccurate, probably because of the spatial resolution of the as- 
sociated images. 

The cluster extension has been estimated considering only 
visual or infrared images, but from an inspection of Spitzer 
3.6,4.5,5.8, and 8.0 /mi images it is easy to see that the cloud 
around the cluster is more extended. Spectroscopy of stellar ob- 
jects, covering not only the central region of Sh2-152 but the 
area marked with the (a) square (see Fig. [TJ, would allow us 
to determine their individual distances and the extension of the 
cluster. 

With the mass of associated molecular clouds and high-mass 
protostellar obje cts (HMPO), the mass of Sh2-152 can be indi- 
rectly estimated. lAo et all d2004l) determined the dust mass of 
the CS, 13 CO, and C ls O cores, mapping their respective lines. 
Assuming a distance of 3.5 kpc, from the column densities of 
the three observed molecular lines, they measured masses of 
1.78 • 10 3 M o (CS core), 9.20 • 10 3 M Q ( 13 CO core), 0.29 • 1O 4 M , 
and 1.20 • 10 4 M o (for two C ls O cores). 

From the CO column den s ity and adopting a distance of 5.3 
kpc dWouterloot et all 1 19931) . iGuan et al.l (120081) determined a 
mass of 3.81 ■ 10 4 M o for the associated CO core, from which 
they observed the emergence of mass outflows. The size of the 
CO cloud deter mined by them was 4. 1 pc, similar to the 4.3 
pc reported by lAo et all d2004l) to the 13 CO core, but larger 
than the other three observed cores (1.8 pc for the CS core 
and 1.7-3.0 pc for the C ls O ones). Using 850/iin polarimetry, 
ICurran et all d2004l) studied the region. For IRAS 22566+5828, 
they found outflows and a high-mass protostellar object (HMPO) 
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Fig. 1. Spitzer 3.6 /mi image. It covers the target (a) and control 
(b) fields for Sh2-152. Red squares are equivalent to the LIRIS 
field of view (4.2 arcminutes length). 



with an estimated total mass (dust) of 3500 M Q . In the case 
of IRAS 22566+5830, they detected a magnetic field aligned 
north-south. In their work, they used 5.0 kpc as the distance to 
the object and 100 as the gas-to-dust ratio. 

As the physical parameters depend on the square of the dis- 
tance, it is crucial to have a reliable determination of the dis- 
tance. We therefore refine the distance and extension of Sh2-152 
using near-infrared photometry for the field and spectra for 13 
stars. This would provide the first distance and extension esti- 
mates using more than one spectrum for candidates of the cluster 
stellar population. Likewise, we plan to obtain the total cluster 
mass by using its stellar population and luminosity function. 

This paper is organized as follows. Section 2 details the ob- 
servations, as well as the related data reduction procedures. The 
analyses of the photometry and spectroscopy (spectral classifi- 
cation) are described in Sect. 3. The cluster parameters are dis- 
cussed in Sect. 4. Finally, Sect. 5 summarizes our conclusions. 



2. Observations 

Our study is based on broad-band infrared imaging (J, H, 
K$) and medium resolution multi-object spectroscopy (H and 
K), acquired with LIRIS, a near-infrared imager/spectrograph 
mounted at the Cassegrain focus of the 4.2 m William Herschel 
Telescope (Roque de Los Muchachos Observatory, La Palma). 

This infrared data was complemented with optical spec- 
troscopy for the central ionizing source of the cluster, obtained 
with FIES at the Nordic Optical Telescope (NOT, Roque de los 
Muchachos Observatory, La Palma). A summary of the observa- 
tions is given in Table Q] 
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Table 1. Summary of imaging and spectroscopic observations. 



RA (J2000) Dec (J2000) 


Filter 


Exp. time 


Seeing 


[km.] [° ' "] 




M 


["] 


Sh2-152 field imaging: 


22 58 37 +58 46 24 


J 


222.30 


0.62 




H 


93.12 


0.69 




K s 


66.12 


0.58 


Multi-object infrared spectroscopy (R ~ 


2500): 




22 58 45 +58 46 50 


H 


2400 


1.04 




K 


2400 


0.98 


Optical spectroscopy (R ~ 46000): 


22 58 40.8 +58 46 58.2 


370-730 nm 1555 


1.50 


Control field imaging: 


22 58 50 +58 38 24 


J 


79.34 


0.62 




H 


79.34 


0.59 




K s 


79.34 


0.65 



2.1. Imaging 

The cluster candidate images were obtained on 2006 July 22 
and the control field was observed on 2009 November 27. In 
both runs, seeing was between 0.58" and 0.69". Observations 
were carried out using LIRIS, an infrared camera equipped with 
a Hawaii 1024 x 1024 HgCdTe array detector, with a field of 
view of 4.2' x 4.2' and a spatial scale of 0.25" pixeV x , using 
J (A c = 1.250 /mi, A/1 = 0.1 60 /mi), H (A c = 1.635 /mi, 
AA = 0.290 //m), and K s (A c = 2.150/mi, AA = 0.320 /art) 
filters. To improve cosmic-rays and bad-pixel rejection, and to 
construct the sky image for the sky subtraction, we observed 
in eight-point dithering mode. Data reduction (bad pixel mask, 
flat correc tion, sky subtract i on, an d alignment) was done with 
FATBOY (Eik enberrvetall 120061) and geometrical distortions 
were corrected with the LIRIS reduction package, LIRISDPQ. 
The final false-colour image for Sh2-152 field is shown in Fig. 

m 

The instrumental photomet ry was done w ith DAOPHOT II, 
ALLSTAR, and ALLFRAME dStetsonl [19941) . and was cleaned 
of non-stellar and poorly measured objects based on two param- 
eters. The first one is the sharpness index, which is the difference 
between the observed width of the object and the PSF width. 
Objects such as hot pixels have a negative sharpness values while 
a "broad" object, such as galaxies or pairs of unresolved stars, 
have a positive sharpness. The sharpness was complemented 
with the estimated standard error in the measured magnitude, <x. 
Using values of sharpness in the range -0.25 < sharp < 0.25, 
and cr < 0.1, we filtered 20% of the detections, all of them with 
Kg over 16. 

The photometric calibration (inc luding a colour term) was 
done using the 2MASS catalogue ( ISkrutskie et all 120061) for 
the three filters. We selected isolated and non-saturated stars, 
with photometry in 2MASS well-distributed over our whole im- 
age and the colour-magnitude diagram. For the stars with Ks 
brighter than 1 1, we adopted the 2MASS photometry, owing to 
their image saturation. 

The astrometric calibration was performed with SKYCAT, 
matching physical and equatorial coordinates for 21 sources, 
well-distributed in our LIRIS image. Typical values of the RMS 

1 http://www.iac.es/project/LIRIS 



of the fitting were smaller than 0.15 arcsec for the three bands, 
which is adequate for the mask design. 

2.2. Infrared spectroscopy 

The infrared spectroscopic observations were obtained on 2008 
August 19 with seeing ~ 1" and using the multi-object spectro- 
scopic mode (MOS) and the HK pseudogrism, with a resolution 
of A/AA ~ 2500. We observed two AO V stars as telluric stan- 
dards: HD 223386 for the H pseudogrism observations and Hip 
107555, for the K one. The mask contained 13 slits of width 0.8 
arcsec and length 6 arcsec and was centred at o^ooo = 22 /! 58'"45 J 
and (52ooo = +58°46'50". In Fig. [2] we mark the stars observed 
spectroscopically in the near-IR. The selection criteria for multi- 
object spectrocopy is detailed in section [XT] 

We observed using an ABBA strategy (a star is located in 
positions A and B, in the slit, and then sequentially changed); 
with this mode, we are able to remove the sky from the spectra 
by simply subtracting spectra at position A from spectra at po- 
sition B, and vice versa. Flat-fielding, tracing, sky subtraction, 
coaddition, and extraction were done using LIRISDR, a package 
developed specifically for LIRIS data that uses the information 
from the mask design files. Thanks to the combination of the 
individual spectra (eight for each band), we discarded cosmic 
rays and hot pixels that could mimic spectral lines. For wave- 
length calibration, argon and xenon lamps were observed, using 
both lamps (continuum subtracted) to calibrate the /T-band spec- 
tra and the argon lamp for //-band calibration. 

Telluric subtra ction was done using XTELLCOR 
(IVacca et all 120031) . an IDL program that, applying a 
high-resolution synthetic model of an A0V star (the spectral 
type of our standard) over the observed telluric standards, 
produces the calibration spectrum with the telluric lines. This 
spectrum was then used to correct our science spectra, with the 
IRAF3 task TELLURIC. 

2.3. Optical spectroscopy 

The optical spectra for the central ioniz ing source of Shl52-2 
(star number 4 from lRusseil et all 120071) was obtained on 2009 
November 12, using the FIES cross-disperser high-resolution 
echelle spectrograph, at the 2.5 m NOT. We used FIES in 
medium resolution mode (R = 46000), with a fibre size of 1.3" 
and the spectral range of 370 to 730 nm covered in a single fixed 
setting. The exposure time for the individual spectra was 25.83 
minutes, and the seeing during the observations was 1.5 arcsec- 
onds. 



3. Results 

Previous near-infrared studies of Sh2-152 were based on the 
2MASS photometry, thus their results were limited to magnitude 
Ks ~ 13. With our new observations, we obtained near-infrared 
photometry with limiting magnitude Ks ~ 20 for the cluster and 
control field. These data were used to select the OB-type star 
candidates from the cluster main-sequence for the spectroscopic 
follow-up. 



2 IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with the 
National Science Foundation. 
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HA (J2D00) 



Fig. 2. LIRIS false-colour image of Sh2-152 (blue=/, green=//, red-Ks); marked with green squares, the stars observed spectro- 
scopically in the near-infrared. In the centre, indicated by a red circle, the central star observed spectroscopically in visual. The 
positions of the masers within the field are indicated by yellow crosses. The (orange) ellipses correspond to the l-<x position uncer- 
tainties for the IRAS 22566+5830 (central) and IRAS 22566+5828 (lower) sources. Black and white colours for central stars are 
used for visibility reasons. 



3.1. Colour-magnitude diagram 

The cluster candidate and control field colour-magnitude dia- 
grams are shown in Fig. [3] In these diagrams we see how the 
cluster's main sequence is scattered following the reddening vec- 
tor due to the differential reddening over the field. This scatter 
also affects part of the local dwarf sequence, which is the prin- 
cipal feature in the control field CMD. Because of the differen- 
tial extinction over the candidate field, cluster stars have redder 
(J - Kg) colours and dimmer K$ magnitudes. 

The selection of OB-type candidates considers this effect by 
selecting stars with (J - Ks) in the range 0.5 - 4.5 and K$<13. 
This first requirement would select stars embedded in the clus- 
ter's natal cloud or their own particular natal cloud, while the 
second requirement comes from the instrumental limitation of 
spectra of high signal-to-noise (SNR). It is necessary to keep 
in mind that telluric variations impose an upper limit on the in- 
frared exposure times. 

For the candidate selection, we also took into account the 
positions of these stars over the field of view. First of all, the 
cluster's boundaries do not appear to be clearly defined in the 
LIRIS near-IR images. Including stars that are not only in the 
central region of the image would help us to derive information 
about the membership of stars separated from the central part 
of Sh2-152. Secondly, the number and positions of slits in the 



mask is limited by spectral overlap in the dispersion axis and the 
spectral features that we wish to observe. 

From the target and control field CMDs, it is also possible 
to estimate roughly the cluster candidate distance. Observing 
both CMDs, we can identify the local dwarf sequence, formed 
by dwarf stars situated in the solar neighbourhood, as a vertical 
strip between (J - Ks ) = 0.4 and {J-Ks) = 1 .0. This strip looks 
sharper for the control field CMD in the whole Ks magnitude 
range but, in the target CMD, this sequence spreads for K$ >15. 
The effect can be understood as the extra reddening caused by 
Sh2-152 over the disc stars located behind the cluster candidate. 
Hence, any star from the local sequence brighter than the magni- 
tude where this sequence widens would be located between the 
cluster candidate and us. By locating this point, we can make a 
first estimation of the distance to Sh2-152. 

A method of obtain ing this estimate is presented in 
iMarin-Franch et al.l d2009t) . Considering 11.2 + 0.75 Gyr as 
the character i stic a ge for stars in the solar neighbourhood 
dBinnev et all l2000h . in their CMDs they plot isochrones with 
this age at different distances and thus under different extinc- 
tions. They then locate the point where target and control field 
local sequences start to differentiate into this set of isochrones. 
The closest isochrone to this point indicates a rough distance es- 
timate for the cluster. 
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Fig. 3. Calibrated colour-magnitude diagrams for Sh2-152 (a) and the control field (b). The black lines represent isochrones of solar 
metallicity and 1 1.2 Gyr, located at three distanc es (3.0, 4.0 and 5.0 kpc) and the grey curve, the positio n of a K2 III star at s evera l 
distances. The isochrone was generate d using thelPietrinferni et al. l d2004l) stellar evolutio nary library, the Caste Hi & Kuruczl (120031) 
bolometric correction library an d the jRieke et aiT dl989|) extinction law w ith R = 3.09 jRieke & Lebofskvi Il985h . For the K2III 
sequence, My was obtained from lCoxl (l2000h and IR intrinsic colours from lDucati et al.l (1200 ll) . 



After including three isochrones in our CMDs for 3.0, 4.0, 
and 5.0 kpc, we can see that differences appear between the 3.0 
and 4.0 kpc isochrones (Fig. [3]). Therefore, a first estimate of the 
distance of Sh2-152, completely independent of the spectropho- 
tometric distance estimate, would be less than 4.0 kpc. 



over 97% of the CMD. The value of the completeness factor ac- 
cording to the position in the CMD is shown in Fig. [4] It can be 
seen in the top diagram that for Ks fainter than 16.5, the com- 
pleteness drops to 0.8, meaning that 80% of the stars would be 
detected at that magnitude. 



3.2. Completeness tests 

To precisely remove the disc stellar population from the cluster's 
luminosity function, and quantify the instrumental sensitivity to 
stellar magnitudes and colours, we carr ied out a completeness 
test, as described in lHidalgo et al.l (12008b and lAparicio & Gallar j 
(fl99l . 

Completeness tests consist of adding artificial stars to the 
images and then, repeating the photometry using the same pa- 
rameters and procedures as for target and control field photom- 
etry. For our tests, we included 480200 artificial stars over 200 
images (for J and Ks filters) with Ks instrumental magnitudes 
between 12 and 22 and colour (J - Ks) between and 6; with 
these limits, we homogeneously covered the whole range of the 
CMDs. To avoid the creation of extra crowding, the stars were 
added following a network, within the full range of coordinates 
into the image, and using a star-to-star distance of two times the 
PSF radius plus one. In this array, star positions are fixed from 
image to image, but J and Ks magnitudes vary. 

Finally, for each artificial star of magnitude Ksj and colour 
(J-Ks)i, we defined a section of K s ,\ +0.4 and ( J - K s )\ +0.4, 
where we counted the number of recovered and injected stars. 
The ratio of these numbers is the completeness index assigned 
to that star. The completeness test reproduces the detection limit 
for the CMD and gives a completeness of more than 0.90 for 



3.3. Near-infrared spectral classification 

We based our near-I R spectral classificati ons for the OB stel- 
lar types on t hose o f lHanson et al. I (fl996h for the /f-band and 
lHanson et al.l (119981) for the //-band, a nd for the later spec- 
tral typ es on those of iMever et al.l dl998l) and I Wallace & Hinklel 
dl997l) . The characteristic lines for each spectral type were com- 
plemented with a visual comparison between our spectra and 
other spe ctral catalogues, with the H- and /T-bands at similar res- 
olutions (llvanov et al.l, 120041 iRanadee t al., 20041 iRanada et all 
120071 lHanson et all l2005h . Table [2] contains the coordinates, 
near-infrared magnitudes, and spectral types of the spectroscopi- 
cally observed stars. The positions in the CMD and CC diagrams 
for the observed stars are presented in Figs.[5]and[6] 

Spectral classification in the /f-band can be problematic. 
The main feature in this band for late B to early F-type stars is 
Brackett y at 2. 166 /mi, but this line is adjusted in the AO V mod- 
elling when constructing the telluric-line correction spectra. Any 
over- or under-correction could therefore add noise to this partic- 
ular line. Complementing the analysis with the //-band spectral 
information overcomes this problem. 

Another consideration stems from the nebular differential 
emission that dominates the central zone of the cluster. This af- 
fects the stellar spectra in two ways. First, the emission can con- 
taminate the spectra, making it hard to distinguish the nature of 
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Fig. 4. Surface formed by the points in colour-magnitude- 
completeness space is shown. The top figure corresponds to 
the Sh2-152 field and the bottom figure, to the control field. 
Completeness ranges from to 1 . 

the lines, whether stellar or nebular. This is noticeable in star #8, 
whose wide Brackett-series emission lines seem to be contam- 
inated with a narrow absorption line. This absorption line is an 
artefact resulting from the nebular contribution included in the 
slit and reversed through the sky subtraction. 

The second problem is that, owing to the inhomogeneous 
nebular distribution, for the same slit we can have sky with an 
added nebular contribution relative to the local sky or even no 
nebular contribution to the star position. This causes an artificial 
subtraction of nebular lines or, in some cases, an artificial emis- 
sion line. This effect is easy to detect through a visual inspection 
of the individual A-B or B-A spectra, detecting lines that pass 
from absorption to emission as we change from an A-B to a B-A 
image. Examples of spectra with this behaviour are numbers 7 
and 10, both located in the central region of the cluster, where 
inhomogeneous nebular emission is dominant. The final near-IR 
spectra are shown in Fig. [7] 

Star #1 is located in a very interesting sector of the field. Its 
surroundings look to be clear in the /-band but for Ks a bright 
arc and a nebulosity arises. This arc even rivals the central zone 
of Sh2-152 in brightness (Spitzer ch-4, 8.0yum). Because of this 
nebulosity, star #1 is the most reddened star in our set but we 
cannot determine whether it is the ionizing star of the IRAS 
22566+5828. The star was classified as Bl V, based on the shape 
and depth of its Bry line and the presence of helium (2.1 1 fim). 




2 4 

(J-Ks) 

Fig. 5. Calibrated colour-magnitude diagram for Sh2-152. The 
red segmented arrow shows Ay = 20 mag (Ak s = 2.04) and the 
main sequence is located at the distance determined in this study 
(3.21 kpc). The continuous blue line in the upper part of the main 
sequence corresponds to the best fit made to project the CMD 
stars to the main sequence, following the reddening vector. The 
position of the spectroscopically observed stars is marked with 
blue numbers. 

The //-band is quite noisy, and the He 1 1.70 fim is contaminated 
by nebular emission. Owing to this noise, produced by its lower 
position in the detector causing the lost of half of the flux, the 
stellar type could easily be two subtypes earlier or later. 

Spectrum #2 is characterized by its strong emission Brackett 
lines. This object is located below the main-sequence line in 
the colour-colour diagrams (Fig. |6). Its near-IR colour excess, 
which is characteristic of young stellar objects, and the dominat- 
ing emission lines over the whole spectra lead us to classify it as 
a YSO. 

Stars #3 and #4 have similar features and their spectral types 
vary in the range B2-3V according to their //-band spectra. 
Star #3 spectra is noisier than star #4, and this explain the 
weaker Brackett series. However for both spectra the series are 
present until Brl5 (1.57 /zm), confirming the similar spectral 
types of both stars. We also detect for both stellar spectra the 
He I 1 .70 jum line. The Zf-band spectrum does not clearly dis- 
play the He I 2.11 fim line, but the Bry line is again consistent 
with a B2-3 V spectral classification. 

For star #5, the spectral type is difficult to define. The spec- 
trum clearly shows the lines Brl2, Bry, and He I at 1.70/itn but 
a weak Brackett series. This is characteristic of late O- to early 
B-dwarfs. Comparing with spectra from the Hanson catalogues, 
we find that the spectral type of this object corresponds to B2 V. 

In the case of stars #6 (H = 13.75 and Ks = 12.66) and 
#7 (H = 14.03 and Ks = 13.24), it was impossible to derive 
a precise spectral type owing to the presence of some emission 
lines, the lack of other strong spectral features, and the low SNR 
of our data for these dim objects. 

In the centre of the cluster, nebular emission is an important 
contributor to the spectra and, in some objects, appears to be in- 
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Fig. 7. Individual spectra for the observed Sh2-152 stars. The spectra are arranged, from top to bottom, into YSOs, early types, late 
types, and those without spectral classification. The final red spectrum corresponds to the telluric correction. 



homogeneous. Star #8, practically in the middle of the field, has 
been resolved thanks to our LIRIS image resolution. The nebu- 
lar emission is overcorrected because of the inhomogeneities in 
the sky emission around the object. The hydrogen lines exhibit 
a wide emission component mixed with shallower absorption. 
This might be due to the overlapping of the nebular (wide) con- 
tribution with stellar (narrow) lines. We also detect He I at 1.7 
and 2.06//m for this object. The position of the star #8 in the 
colour-colour diagram, below the main-sequence line, indicates 
that it is a YSO-like star #2. Because it is less reddened, it might 
have cleared its surroundings. 

A similar effect could occur for star #9; its spectra looks fea- 
tureless in absorption, probably because its lines are filled by the 
nebular emission lines. The classification of the spectrum, dom- 
inated by Brackett and He I emission lines, was impossible. 

Star #10 has a well-defined Brackett series, with no indica- 
tion of either helium or carbon lines. A comparison with star 
HR5685 from Ranade's library leads us to classify this object as 
B5V. 

Stars #11 and #12 are late-type stars, even if the CO bands 
do not clearly appear in these spectra. A close similarity can 
be found between star #1 1 and Ranade's HR5019 (spectral type 



G6.5 V) and between star #12 and Meyer's HR4496 (spectral 
type G8 V). Both appear to be stars in the Galactic disc. 

Finally, star #13 has a remarkable position in the colour- 
magnitude diagram, being the brightest star in our selection. 
However the weak Brackett and the strongest CO heads in our 
candidates indicate a late-ty pe for this star. The spectra looks 
similar to those of HR5888 dlvanov et all |2004l) and HR7328 
dMeveretall[l998t) . We classify star #13 as G8-G9III. 

3.4. Spectral classification (optical) 

Our final optical spectrum for the central star of Sh2-152, is 
shown in F ig. [8] This star was previo usly classified as Q9 V by 
ICrampton et aU (fl978l) and 08.5 V bv lRusseil et all d2007b . 

Our spectral clas s ificati on is based on the study of 
IWalborn & Fitzpatrickl d!990h . On first inspection, the spectral 
type for the central star of Sh2-152 appears to be O9-B0 V star, 
because of its ratio of He I 4471 A to Hell 4542 A. However, 
there are a few spectral features in the optical spectra that de- 
serve some review. 

First of all, the He II 4686 A line is similar in depth to the 
He I 4471 A line and deeper than the He I 4388 A one. The 
ratio Hell 4686 A-Hel 4471 A looks slightly larger than for 
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Table 2. Spectroscopically observed stars. 
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12.127 
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+58 46 59.56 


14.395 


13.343 
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10 


22 58 40.422 


+58 47 06.06 
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12.863 
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B5 V 


1 1 


22 58 42.920 


+58 47 26.67 


13.572 


12.614 


12.010 


G6.5 V 


12 


22 58 48.062 


+58 48 00.22 
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13.923 
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G8 V 


13 
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+58 48 36.34 


11.203 


10.491 


10.344 


G8-9 III 
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Fig. 6. Calibrated colour-colour diagrams for Sh2-152. The blue 
numbers corresponds to the spectroscopically observed stars. 
Red arrows in both diagrams show the reddening vector. 

stars between 09 V and B0 V. An He II 4686 A spectral line 
stron ger than that cen tred on He I 447 1 A would indicate a Vz 
type ( IWalbornl 120091) and might indicate that there is ongoing 
massive-star formation in the cluster. 

Secondly, the CNO complex at 4630-4650 A looks slightly 
fainter than expected for those spectral types. Although this 
might be caused by the continuum rectification, this would not 
seem to be sufficient to justify the difference. 

The third interesting feature is the absence of a Si III (4553- 
4574 A) detection. These lines appear for spectral types later 
than 09.5 V, hence this detection indicates a spectral type of 
09 V or earlier. However, for types earlier than the aforemen- 
tioned limit, the depths of the Hell 4542 A and Hell 4686 
A lines do not fit. 

Therefore, the star was assigned a spectral type of O9-B0 V, 
with spectral features that are not fully consistent with any of 
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Fig. 8. Optical spectrum for the central star of Sh2-152 
(Russeil's star number 4). Vertical segmented lines (grey) mark 
the spectral lines used in the spectral classification. Overlapping 
zones of the spectral orders are marked with red vertical lines; 
the low SNR causes a poor combination of the spectral orders, 
which could mimic a spectral line. Spectra resolution was de- 
creased to R — 8000, to help us to discern the lines. 



the subtypes. While the SNR of the spectrum is relatively mod- 
est (the SNR per pixel is 20, but at R = 46000, which allows us 
to degrade our spectrum without loss of information), we can- 
not exclude the presence of a binary system, a composite spec- 
trum, or a slightly peculiar spectrum. The possibility of a Vz 
type would also enhance the peculiarity of this star, but data of 
higher S NR is mandatory to confirm this point. 



4. Discussion 

4.7. Distance and size determinations 

Using the estimated spectral types, we derived the individual dis- 
tances tocluster stars. Assuming the absolute visual magnitude s 
of ICoxl (120001) . the intrinsic infrared colours of lTokunagald2000l) . 
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Table 3. Reddening and distances of early-type stars. 



Spectral type 


ID 


Ak s 


Distance 






[mag] 
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Bl V 


1 


2.70 


3.17 


B2 V 


5 


0.70 


3.76 


B2-3 V 


3 


0.76 


3.61 


B2-3 V 


4 


0.62 


2.70 


B5 V 


10 


0.95 


2.83 



and the lCardelli et al.1 d 19891) extinction law with R = 3.10, the 
extinction for the Ks band can be expressed as 
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Fig. 9. Radial distribution of stars for Sh2-152 (solid black 
curve) and the control field (red dashed curve). 



Across the field of view, we can find several situations of extinc- 
tion ranging from 0.62 to 2.70 mag for Ak s , or Ay between 5.44 
and 23.68 mag. In Table[3] we present the values of Ak s and the 
individual distances of the early-type stars used in the distance 
estimation. 

We excluded Russeil's star number 4 (central star of Sh2- 
152) because it appears saturated in our images, and 2MASS 
images of this object are unable to resolve the central zone of 
the cluster. The 2MASS magnitudes for the central star would 
be contaminated and we prefer to exclude this object from the 
distance estimation. 

With the five early-type stars, we obtained an average dis- 
tance to the cluster of 3.21 + 0.21 kpc from the i ndividual cal- 
culatio ns. This value is larger than that given by iRusseil et alJ 
d2007l) - D — 2.39 kpc - and slightly smaller t han the distance 
of 3.5 kpc estimated by ICrampton et all dl978l) . Both estimates 
were obtained from spectral observations of the central ioniz- 
ing star of the cluster (our O9-B0 V star). In our case, the es- 
timate was derived from individual distance estimates that are 
self-consistent and cover a wider area of the cluster. 

From the Spitzer image of the field and our observation of 
a nearby control field using LIRIS, we learned that the angu- 
lar size of the cluster is comparable to the LIRIS field of view 
(4.2' X 4.2'). The cluster was previously considered as just the 
surroundings of Sh2-152, underestimating the size of the object, 
but from our CMD, we realize that stars near the edge of our im- 
age are part of the cluster. Hence, considering at least a cluster 
size of 4.2' (the LIRIS field of view) and the average distance of 
3.21 kpc, the diameter of the cluster would be > 4.0 pc. 

A test that helps us to clarify the extension of the cluster is to 
count the number of stars for different concentric rings covering 
the whole field of view. Taking as centre the coordinates of IRAS 
22566+5830, we counted the sources in concentric rings of 25 
arcsec in thickness. In Fig.|9j we can see that there is an evident 
change in the slope of the counts in the central section of the 
image (around 250 pixels, i.e. 60 arcsec), but even as the number 
of stars detected continues to decrease it remains larger than the 
number of stars detected for the control field. In addition, the 
existence of stars with the same spectrophotometric distances as 
objects within the central area of the cluster allows us to confirm 
that the size of the cluster is at least comparable to our field of 
view, and that the subtraction of stellar counts using the same 
field is unreliable. 



4.2. Cluster mass and age estimations 

To derive the cluster mass, we integrated the stellar mass func- 
tion of Sh2-152, previously corrected for the field stellar contri- 
bution, after accounting for the effect of completeness on the star 
counts. 

The first step in obtaining the mass function is to project 
every star following the reddening vector, to the d warf s tar se- 
quence defined by the magnitudes and colours of ICoxl (I2000h 
and located at our estimated distance. This sequence can be ex- 
pressed analytically by a second degree function or a set of first 
degree functions. We found that fitting the dwarf sequence us- 
ing two lines, one from 09 V to AO V and the second from AO V 
to F0 V (blue solid lines in Fig. [5j was more effective than fit- 
ting a second-degree function. The cut in F0 V arises from our 
detection limit; at later spectral types (fainter magnitudes), the 
histograms start to decline towards lower masses. 

Once the stars in the Sh2-152 and the control field CMDs 
were projected following the reddening vector to the line-fitted 
main sequence, we produced the histograms including the spe- 
cific weight for each star, according to its completeness factor 
determined by its position in the CMD (J - Ks , Ks). The correc- 
tion for the incompleteness consists simply of multiplying the 
number of stars by the inverse of the completeness factor. 

The histograms bins were 0.8 in Ks, equivalent to twice the 
size of the completeness bins. To convert a Ks magni tude to a 
stellar mass, we used the values given bv lCoxl(l2000l) . For in- 
between magnitudes catalogued values, we interpolated between 
the two closest values. 

After subtracting both histograms (i.e. subtracting the con- 
trol field histograms from the Sh2-152 field), we derived the 
cluster mass function (Fig. [TOt . Deriving the best-fit linear re- 
lation (F = -1 .62) in the log(N)-log(M) plane and then integrat- 
ing from 0.20 to 20 M , we obtained a total mass for Sh2-152 of 
(2.45 + 0.79) ■ 10 3 M Q . 

To compare this result with the literature, we had to care- 
fully consider the nature of both quantities. Our value comes 
from evaluating the st ellar con t ent of Sh2-152, in sharp con- 
trast to studi e s such as lAo et all d2004l) . iGuan et all (l2008h . and 
ICurran et al.l d2004h . who derived the virial mass, i.e. the stellar 
mass necessary to ensure that the molecular core is gravitation- 
ally bound, from molecular radial velocities. 
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Fig. 10. Mass function for Sh2-152. The points indicate the cen- 
tral position of the mass range indicated above, equivalent to K$ 
magnitude bins of 0.8. The continuous red line shows the best 
linear fit to the data. 



iGuan et al. | J2008I) calculated the mass (dust) for the asso- 
ciated CO core. Their value for the total mass (dust) is 3.81 ■ 
10 4 M o , an order of magnitude higher than our resul t. If we use 
our di stance estimate instead of the value used by IGuan et al.l 
(120081) (5.3 kpc), the total mass (dust) for the CO core drops to 
1.22 • 10 4 M Q and its virial mass is also reduced by a factor of 
0.32, down to 3.04 ■ 10 3 M p. 

The study of lAo et alJ (12004 assumes a distance that is more 
similar to our estimate. With d = 3.5 kpc, they obtained four 
total mass (dust) estimates: 1.78 • 1O 3 M for the CS core, 9.20 • 
10 3 M o for 13 CO core, and (0.29 - 1.20) ■ 10 4 M Q for the two 
C I8 cores. With the exception of the virial mass derived from 
the 13 CO estimate (4.90 • 10 3 M o ) , their assoc iated virial masses 
are lower than those obtained bv IGuan et al J 1.78 - 10 3 M Q (for 
CS core) and (0.67 - 1.40) • 10 3 M G (for the C ls O cores). 

With the stellar mass and dust mass estimates, we can eval- 
uate the cluster's star formation efficiency (SFE). The star for- 
mation efficiency is the ratio of the cluster's stellar mass to the 
gas-stellar mass. To derive the total gas mass associated with 
the dust total mass, we used the gas-to-dust ratio. According to 
ICurran et al ] (12004 . this ratio varies from 45:1 to 100:1, so we 
prefer to use both values in the estimate; 

Using the dust mass measured by ICurran etal ] (12004 . pro- 
jected to our estimated distance instead of the 5.0 kpc used by 
them, we obtained the associated gas mass, and estimated the 
SFE to be b etween 1.7% and 3. 7%. This value is similar to that 
reported by iHunter et al.l d 19901) (i.e. 6%), assuming a distance 
of 4.0 kpc and estimating the stellar total mass with a normaliza- 
tion of the IMF to the upper stellar range, which in their case is 
18-25M Q . 

We can neither precisely determine the age o f the cluster 
nor c onfirm the idea of sequential stellar evolution (IChen et all 
I2009h . but we can constrain the cluster's age based on the physi- 
cal characteristics of its stellar population. First of all, our earli- 
est main-sequence star is of type 09 V. Since it has not evolved 



away from the main sequence, the cluster lifetime is shorter than 
the time spent by this star on the main sequence. Assum ing that 
18 M Q is the mass of an 09 V star dMartins et aUl2005l) and that 
the cluster population was formed at the same time, we can only 
limit the cluster's age to be younger than 9.4 Myr, using the stel- 
lar mo dels for a 20 M Q solar metallicity star of ISchaller et al.l 
(Il992h . This limit is obviously conservative and cannot be con- 
sidered an age estimate. 



4.3. Ionizing sources 

We used our spectral type classification together with data col- 
lected from the literature to estimate whether the total ioniz- 
ing photons emitted by the stars can explain the radio obser- 
vations of Sh2-152. To determine the number of ionizing pho- 
tons (log Qo) for O-type stars, we used lMartins et al. (120051). For 
B0 V s tars, we used the values given bv lSimon-Dfaz & Stasihskal 
(120081) and for the B0.5-B 1 V spectral types we took the radii 
from the binary CW Cep dBozkurt & DegirmenciL 20071) . the 
spectral type-r e fj relation from iLefever et al' I (120101). and the 
T e s - log qo from the B-star grid from lLanz & Hubenvl d2007l) . 

To derive the number of Lyman continuum ph otons from 
the radio density flux, Nl y c, we used Eq. A. 5 of iHunt et al.l 
d2004 . assuming an electron temperature of T — 10 4 K and 
«(He + ) «: n(H + ). Since this formula is for radiation-bounded 
ionized H II regions, it is not guaranteed that it is applicable to 
IRAS 22566+5828, although an estimation may be illustrative. 

We divide the ionizing stars into two groups: the first con- 
tains stars in the central region (r ~ 30"), including IRAS 
22566+5830, star #10 and Russeil's star number 4. The earliest 
star in this region is Russeil's number 4 (O9-B0V, this work), 
with a total number of ionizing photons between log Qo = 48.06 
and 47.29. The number of Lyman continuum photons derived 
from t he radio observations for IRAS 22566+5830 is log N^yc = 
48.02 (ICondon et al.L [l998l) . which is most consistent with an 
09 V spectral classification for the central star of Sh2-152. 

The second region, formed by stars #1, #3, #4, and #5, is lo- 
cated close to IRAS 22566+5828 source and its earliest star is 
star #1. Within our sample, star #1 would therefore be the main 
ionizing source associated with IRAS 22566+5828. However, 
the value of log Qo - 45.66 for this star is significantly lower 
than for the Lyman continuum ph otons, logA ^ Lvc = 48.35, de- 
rived from the radio observations dHariu et al.Lfl9 98) indicating 
that the star #1 does not produce enough photons to explain the 
radio observations. The K$ image shows a deeply embedded star 
in the centre of the red nebulosity close to of IRAS 22566+5828, 
which appears to be the 07.5-8 V star required as the single ion- 
izing source of this region. Future K spectroscopy for this star 
would help us to confirm this result. 



5. Conclusions 

In the framework of our MASGOMAS project, we have com- 
pleted a spectrophotometric study of the stellar content of the 
compact HII region Sh2-152 using LIRIS at the WHT and the 
FIES cross-disperser echelle spectrograph at the NOT. 

Our new near-infrared data have allowed us to resolve the 
central part of the cluster, to select OB candidates for H and 
K spectroscopy, and estimate their individual distances. The 
mean of the individual distances to five B-type stars, 3.21 ± 
0.21 kpc, pl aces the cluster sl ightly c loser than the estimates 
obtain e d bv ICrampton et"aTI d!978). IWouterloot & Walmslevl 
(ll986l) . IWouterloot et all dl988h . and faariu et all d 19931) but def- 
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initely further away than estimated by iRusseiTe t al. (2007|). In 
any case, our estimated distance is in closer agreement with the 
distance estimates based on the central ionizing star (star num- 
ber 4 in lRussei l et al. 2007) than with those derived from 12 CO 
radial velocities. 

We also obtained the cluster mass function, corrected for 
the Galactic disc stellar contribution (using the observed con- 
trol field) and for completeness effects (through artificial star 
injection and photometry). Integrating this function, we cal- 
culated a lower limit to the cluster total mass. The value of 
(2.45 + 0.79) • 1O 3 M agrees with the mass values from the liter- 
ature, all of them being virial masses necessary to reproduce the 
observed molecular velocities (for both CO and CS cores). 

The central O9-B0 V star remaining on the main sequence 
allows us to limit the cluster's age to be younger than 9.4 Myr. 
The optical spectra for this central star displays several interest- 
ing features. First, the strenght of its helium lines (He II 4686 
A similar to that of He I 4471 A and stronger than He I 4388 A) 
might be indicative of a Vz spectral type, although the relation 
between the He II 4686 A line and He II 4542 A line appears to 
exclude this. Secondly, the spectra detect a fainter CNO complex 
at 4630-4650 A for the derived spectral type. And third, we have 
found the absence of Si III (4553 A and 4574 A). Either a binary 
system, a composite spectrum, or a slightly peculiar spectrum 
are possible explanations of the last two features. Future obser- 
vations of this object with higher SNR will help us to shed light 
on the real nature of the central, most massive star of Sh2-152. 

For two regions of Sh2-152, we have compared the number 
of ionizing photons from the classified OB-stars with the number 
of Lyman continuum photons derived from radio observations. 
For the central region containing both IRAS 22566+5830 and 
the O9-B0 V central star (Russeil's star number 4), both of these 
numbers are consistent with there being a single ionizing source. 

For the second region containing IRAS 22566+5828, the 
number of Lyman continuum photons cannot be explained by 
the earliest star in our sample, and the derived value of log N^yc 
would imply the existence of an 07.5-8V star. This source 
could be a highly embedded star, without spectral classification, 
at the centre of the bright red arc close to star #1. We expect fu- 
ture spectral observations of this star to help us verify this find- 
ing. 
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